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1. Intro 



The dawn of the first astrophysical objects and the eventual reionization of the Universe is 
of fundamental importance, offering insight into many astrophysical processes. Also, since ob- 
servational surveys only pick out the rarest, most luminous objects, cosmic radiation fields and 
their effects on the intergalactic medium (IGM) may our only way of studying the majority of our 
galactic ancestors. Recent theoretical and observational advances are starting to shed light on a 
complex and extended reionization epoch. Here I sketch-out three likely stages of this process. I 
briefly discuss the potential of the redshifted 21cm signal to probe not only reionization, but also 
the pre-reionization epochs. However, interpreting the upcoming observations requires efficient 
modeling tools, such as the recently developed 21cmFAST. 



2. Once upon a time... the story of cosmic dawn and reionization 



2.1 Early stages 

Hierarchal structure formation implies that the first astrophysical objects are likely hosted by 
so-called minihalos (with virial temperatures T^^ < lO'* K). Cooling via atomic hydrogen is ineffi- 
cient at such low temperatures, so baryons condense into these halos via the molecular hydrogen 
(H2) cooling channel. Without metals to aid in cooling and fragmentation, the first generations of 
stars (so-called PopIII stars) are likely massive (with masses of ~ lOOM©; ||, |8^), yet short- 



lived (with lifetimes of a few Myr; |62, 53|). 

PopIII stars are likely to have very different properties from "normal" PopII stars, with harder 
spectra, a factor of ten times more ionizing photons per baiyon [|6^], and various exotic, mass- 
dependent fates (e.g., [30]). Accurately determining their properties, formation efficiency and 
initial mass functions (IMFs) requires detailed numerical simulations and costly explorations of 



parameter space (e.g., []74p). 

These prehistoric giants have complex interactions with their surroundings, since star forma- 
tion in such small-mass halos is susceptible to feedback mechanisms. Modeling these requires 
introducing parameterized prescriptions for stellar properties and extending simulations to larger, 
~ 1 Mpc scales. Feedback processes can come in three forms: (1) mechanical: supemovae (SNe) 
can blow out gas from the surrounding shallow potential well, delaying local star formation until 
gas is re-accreted (e.g., [^, U^); (2) chemical: evolution in the IMF depends on metal enrich- 



ment [ p9[ pq , |64| ], which is very inhomogeneous making it necessary to model the enrichment on 
large scales (e.g., [fTl]]); and (3) radiative: radiation in various bands can either promote (positive 
feedback) or suppress (negative feedback) future star formation. 

Because it acts on very large scales and involves a large dynamic range in energy, it is likely 
that radiative feedback is the most pivotal and complex of the feedback mechanisms. Positive 
radiative feedback can result when the enhanced free-electron fraction from ionizing photons (e.g.. 



\p%) or hydrodynamical shocks [ ]65[ ] catalyzes the formation of H2, thereby enhancing the H2 
cooling channel. Negative feedback can result from heating by ionizing radiation which can photo- 



evaporate gas in low-mass halos (e.g., [11 1). Also, an active background of Lyman-Werner (LW) 
radiation (with photon energies in the 11.1 8-13. 6eV range) can dissociate H2, thus decreasing the 



gas's cooling capabilities (e.g., [29]). 
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The first astrophysical objects can emit radiation in various bands. X-rays, generated from 
accretion onto remnant black holes or massive X-ray binaries, have long mean free paths, but their 
feedback effects seem to be mild [|^, However, they can be responsible for early IGM heating 
and some degree of reionization (e.g., ^ 51 1). Feedback resulting from a LW background is 
by definition negative, and can delay star formation in low mass halos, depending on the strength 



of the background ||38|, ^ Also depending on its strength, a transient (due to the short 
lifetimes of PopIII stars) ionizing UV background (UVB) can have various feedback regimes [58 



|47| , p8| , pq , |75| ], especially in conjunction with the negative feedback effects of a LW background 
[ p7| , p^ , 76]. Despite this complexity, it is likely that the final suppression of star formation inside 
minihalos was regulated by the build-up of persistent UV or LW backgrounds before the bulk 



of reionization [26, 27]. However, the complexity of radiative feedback and/or an early X-ray 



background could precipitate an extended early reionization epoch, which could allow a significant 
amount of reionization to occur at late times (z '--'6), and still be consistent with constraints from 



the Wilkinson Microwave Anisotropy Probe (WMAP) [36] 



2.2 Middle stages 

After the initial, "fireworks" stage of PopIII-driven reionization, longer-lived stars begin to 
establish persistent UVBs. Reionization proceeds in an "inside-out" fashion on large-scales, with 



biased sources driving HII regions to tens of Mpc in size (e.g., [21, 19, 85]). Larger, atomically- 
cooled halos (Jyir > lO'^ K) start dominating the ionization budget. These are likely too faint to 



image directly with current or upcoming instruments (e.g., [61, 81]), although we can get close 
with gravitational lensing surveys (e.g., []68|, 53, 80]). Likewise, because of the inside-out nature 



of reionization, most astrophysical objects during the epoch of reionization find themselves inside 
large HII regions, largely unaffected by the neutral IGM even when it still occupies a substantial 
volume filling fraction. Thus most of the current observational probes of reionization (derived 
from quasar proximity regions, damping wings in quasar and gamma ray burst (ORB) spectra, Lya 
emitter (LAE) number density and clustering properties; see below) are most effective back when 
the Universe is predominantly neutral. 

Aside from these observational challenges, there are substantial theoretical hurdles to over- 
come when studying the middle stages of reionization. Ionization structures during this epoch can 
span volumes over four magnitudes larger than the first-stars simulations discussed in the previous 
section (see Fig. 1). To statistically capture the middle stages of reionization, simulation boxes 
need to be hundreds of Mpc in length, yet still be able to account for ionizing radiation from small- 
mass halos. The required dynamic range is daunting, and state-of-the-art simulations are forced to 
ignore minihalos and photon sinks, while also using approximate prescriptions to assign ionizing 
luminosities to source halos (see the recent review in [^]). More approximate techniques have re- 
cently been developed in order to overcome the challenge of such a large range of scales, including 
sub-grid ^ Q] and semi-numerical [|8|, |0[ ||, ||, ^ |lj models. 

Additionally, our poor understanding of the early Universe means that large explorations of 
parameter space are required for any robust conclusion. Analytic models have been very useful in 
this respect, predicting that HII morphology is not very sensitive to redshift, when normalized to 
the same mean neutral fraction, xhi, and also that sources are more important that ionizing sinks 



in this regime [21, 1%. These predictions were confirmed by the first numerical parameter studies 
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Figure 1: Scales of some early Universe simulations, with box sizes indicated in comoving units. Slices correspond to 
metallicity maps from ^1%, and ionization structure from and (left to right). The left and middle panels roughly 
correspond to the current maximum box sizes of state-of-the-art numerical simulations of the first stars and the middle 
stages of reionization, respectively. This figure highlights the daunting dynamic range of cosmological processes such 
as reionization. In principle, one should continue "zooming-in" down to stellar scales (e.g., | |6^ , po||). 



in [44]. Another parameter study, [49], combined ID hydrodynamical collapse simulations with 
large-scale 3D semi-numerical simulations of the UVB during reionization, and concluded that 
radiative feedback does not play a major role during this epoch. Results of such parameter studies 
allow us to simpUfy various theoretical problems. 



2.3 Final stages 

As reionization progresses, HII regions grow to be larger than the mean free path of ionizing 
photons through them, requiring even larger simulation boxes (see Fig. 1). In this final regime, 
sinks of ionizing photons, such as Lyman limit systems (LLSs) could regulate the progress of 
reionization [18], and its morphological structure ijl^,^. These absorption systems are likely found 
in the filaments of the recently-ionized IGM, where the UVB is the weakest [^]. The final "overlap" 
stages of reionization could be delayed until either absorbers get photo-evaporated (e.g. [^]), or 
the remaining neutral islands get ionized from the inside. The time-scales for these processes are of 
order Az 1. This, combined with the large cosmic variance during this epoch (e.g., [^]), suggests 
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that the overlap stage is much more gradual than was initially predicted by small-box simulations 
(e.g., [^5|]). It is also unlikely that these final stages leave a strong observational imprint, such as a 
sharp rise in the UVB [[l8|]. Despite popular belief, current observations cannot rule-out a very late 
overlap, even at z '--^ 5 [37, 46, 42|. 



2.4 Observational constraints 

Several model-dependent constraints on f hi at z 6 have been derived from various astrophys- 
ical probes such as: (1) the size of the proximity zone around quasars ([79, 6|j, but see [ |5^ , 0, 0]); 
(2) a claimed detection of damping wing absorption from neutral IGM in quasar spectra ([^4}, 



but see []52[]); (3) the «o«-detection of intergalactic damping wing absorption in a GRB spectrum 
([^], but see [||]); (4) the number density and clustering of Lya emitters ([^, |8|, |2|, ||, ||], 
but see | [lO| , |53| , [3T|]); (5) the length of dark gaps in quasar spectra ([23], but see §3.3 in [46]). 



Additionally, a direct, nearly model-independent upper limit on xhi was recently derived from the 



dark covering fraction in the Lya and LyjS forests [42]. 

However, constraining the ionization state of the high-redshift Universe is very difficult with 
current observations, since they are not direct probes of HII morphology. For example, the Lyman 
forests begin to saturate at z > 5, becoming difficult to interpret (e.g., [ p7| , p^ , |, 4^). Similarly, 
the dominant population of ionizing galaxies is ~ two orders of magnitude fainter than current de- 
tection limits; and in any case, galaxies offer a very biased tracer of the cosmic radiation fields. 
Luckily, there is a plethora of upcoming observations which should help advance our understand- 
ing: 21cm tomography, high-redshift IR spectra, wide-field LAE surveys, and improved E-mode 
CMB polarization power spectra. Of these, arguably the 21cm line from neutral hydrogen provides 
the greatest potential for directly studying reionization in the near term. 



3. The 21cm frontier and llcmFAST 



The cosmological 21cm signal uses the CMB as a back-light. The offset of the 2 1cm brightness 
temperature from the CMB temperature, Ty, along a line of sight (LOS) at observed frequency v, 
can be written as (c.f. [pO|]): 



57),(v)=27xHi 1 



(l + 5n 



H 



1+z 0.15 



dvr/dr + H 7 V 10 ^Mh^ 



1/2 



0.023 



mK, 



(3.1) 



where xhi is the hydrogen neutral fraction, Ts is the gas spin temperature, is the optical depth at 
the 21cm frequency Vq, 5ni(x,z) = p/p — 1 is the evolved (Eulerian) density contrast, H{z) is the 
Hubble parameter, dvr/dr is the comoving gradient of the line of sight component of the comoving 
velocity, and all quantities are evaluated at redshift z = Vq/ V — 1 . The above equation also assumes 
that dvr/dr <^ H, which is generally true for the pertinent redshifts and scales. 

The redshifted 21cm line carries with it much information, containing both astrophysical (xhi 
and Ts) and cosmological (last four terms in eq. 3.1) terms (see [1^, 17]). However, in order to make 
sense of the upcoming observations, one must disentangle the various components. Considering 
the scales and uncertainties involved, this is a daunting task. 

Recently "semi-numerical" tools have proven invaluable in overcoming these obstacles. These 
simulations use more approximate physics than numerical simulations, but are much faster. In 



5 



How the first generations of luminous baryons established the X-ray and UV backgrounds Andrei Mesinger 



particular, the recent launch of 21cmFAST | |5T| ] (http://www.astro.princeton.edu/~mesinger/Sim) 
provided the astronomical community a publicly-available tool specialized for fast, large-scale sim- 
ulations of the 3D cosmological 21cm signal. For example, a 256^ realization of reionization with 
21cmFAST is obtained in just a few minutes on a single CPU, with every step agreeing with cos- 
mological hydrodynamic simulations into the quasi-non-linear regime ||^,|8|, Such speeds 



facilitate rapid explorations of the vast parameter space, necessary for interpreting upcoming ob- 
servations. 

Most studies thus far have focused on the reionization term, xm, making various simplifying 



assumptions in eq. 3.1. The 21cm signal during reionization can indeed tell us much about the 
relevant physical processes and intergalactic radiation fields. However, the 21cm signal also offers 
us a precious glimpse into even earlier epochs, via the spin temperature, Ts, term. 



The spin temperature can be written as (e.g. [pO|]): 

(3.2) 



1 ~l~ ~l~ -"-c 

where Tk is the kinetic temperature of the gas, and Ta is the color temperature, which is closely 



coupled to the kinetic gas temperature, Ta ~ T^ [|14|]. From this equation, one can see that the 
spin temperature interpolates between the CMB and the gas temperature. It does this through two 
coupling coefficients: (1) the coUisional coupling coefficient, Xc, which requires high densities and 



is effective in the IGM atz > 40; and (2) the Wouthuysen-Field ([78, 13|; WF) coupling coefficient 
which uses the Lya radiation background and is effective soon after the first sources ignite. If either 
coefficients is high, the spin temperature approaches the kinetic temperature of the gas. Otherwise, 
the spin temperature approaches the CMB temperature and there is no signal. 

The CMB temperature decreases as oc (1 + z). After decoupling from the CMB, the gas kinetic 
temperature decreases as oc (i -i-^)^, before being heated, most likely with X-rays from the first 



sources (see Q15| ] and the discussion therein). Likely before the bulk of reionization, the gas is 
heated to temperatures far beyond the CMB, Tk ^ Ty, and the spin temperature term in eq. 3J no 
longer contributes to the signal. When and where Ts < Ty, the 21cm signal is seen in absorption, 
and when/where Ts > Ty, the signal is seen in emission. 

Therefore, the 2 1cm signal and its inhomogeneity could be a powerful probe, not only of reion- 
ization, but also of the Lya background, IGM heating, and cosmology during the coUisionally cou- 
pled regime before the first sources ignite. Another probe of cosmology would be available at lower 
redshifts (making it easier to detect) if heating completes before reionization. A movie of 21cm 
evolution from z = 250, computed with 21cmFAST, is available at http://www.astro.princeton.edu/ 
~mesinger/2 lcm_Movie.html, together with an explanation of the various interesting epochs. 



4. Conclusions 



• Milestones such as reionization are likely the only practical way of studying the primordial 
zoo of astrophysical objects in the near future. 

• Reionization is likely extended, going through various stages. 

• The cosmological 21cm signal is very rich in information, containing both cosmological and 
astrophysical components. 
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• The range of scales and unknown parameter space is enormous. 

• We need efficient modeUng tools, such as 21cmFAST, to make sense of the upcoming obser- 
vations. 

• We are hving in exciting times ! 
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